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First results for three-loop deep-inelastic structure functions in QCD 
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As a first step towards the complete calculation of deep-inelastic scattering at third order of massless pertur¬ 
bative QCD, we have computed the fermionic {nf) contributions to the flavour non-singlet structure functions in 
unpolarized electromagnetic scattering. We briefly discuss the approach chosen for this calculation, the problems 
we encountered and the status of the project. We show the results for the corresponding anomalous dimension in 
both Mellin-A and Bjorken-a:. Together with the m part of A 3 , our calculation facilitates the complete determi- 
nation of the threshold-resummation coefficients B 2 and D 2 . The latter quantity vanishes in the MS scheme. 


1. Introduction 

The computation of the three-loop contribu¬ 
tions to the anomalous dimensions is needed 
to complete the next-to-next-to-leading order 
(NNLO) calculations for deep-inelastic scattering 
(DIS). These in their turn are required for the 
NNLO determination of the parton distribution 
functions (PDF’s) that describe the quark and 
gluon contents of the proton. The accuracy of the 
experiments [ Q has become such that this order 
of perturbation theory is needed to match it. In 
turn, the PDF’s at NNLO accuracy are required 
for precise QGD predictions for the future exper¬ 
iments at the LHG. 

The method we use for obtaining the three-loop 
structure functions is based on the calculation of 
all Mellin moments as a function of the Mellin 
moment number N. Once either all even or all 
odd moments are known, we can obtain the re¬ 
sults in Bjorken-x space by means of an inverse 
Mellin transformation. This method has been 
employed for structure functions since the early 
days of QCD [ ||, ||. Only the coefficient func¬ 
tions are the two-loop level were first obtained by 
a different approach [ i I I, 0- However, the 
fV-space method was later used to confirm those 
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results [||]. 

Until the calculation of the two-loop anoma¬ 
lous dimensions [00 this method could still 
be applied in a rather direct way. The two-loop 
computation of required more finesse, and 

recursion relations [ , also called difference 

equations, entered the scene. The various mix¬ 
ings of the flavour contributions were solved in the 
calculation of a number of fixed moments at the 
three-loop level [ 00 H- This method turned 
out to be so powerful that not only the anomalous 
dimensions could be obtained, but also the three- 
loop coefficient functions as the beginning of the 
NNNLO computation. Hence it seemed best to 
continue along this path and obtain all moments. 
This avoids the difficult problem of operator mix¬ 
ings and gets us the whole three-loop order in one 

go- 

2. Eight Problems 

In this section we discuss eight of the main 
problems that were encountered. Basically these 
are problems that occur in any high-order cal¬ 
culation. Their solution however depends of the 
particular calculation. 

Problem 1 

For the current calculation the mathematics of 
the answer had to be understood better. Hence 
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first the properties of harmonic sum s [ |l0|, 16, 
00 and harmonic polylogarithms [ 
had to be studied. In addition procedures were 
obtained to go from one to the other by means of 
an inverse Mellin transform [ . All this was 

programmed in FORM [ The answer of our 
calculations will comprise harmonic sums of up 
to weight 6 for the Mellin moments. In Bjorken- 
X space we will obtain harmonic polylogarithms 
with weights up to 5. 


Problem 2 

Next we need a scheme in which all integrals are 
reduced to a set of master integrals. To find 
these reduction algorithms in a way that they 
can handle any integral of the necessary topolo¬ 
gies is much work. We have to deal with 0(70) 
(sub) topologies. For some the algorithms are 
trivial, but for the most difhcult cases it may 
take a few months to derive a good algorithm. 
In general the equations are based on integra¬ 
tion by parts [ scaling equations, 

form-factor analysis [|27[| and some equations that 
fall in a special category because they involve 
a careful study of the parton-momentum limit 
P-P^O. 


Problem 3 

The master integrals can usually be determined 
by difference equations. For the determination of 
the rij part of the non-singlet structure functions 
we needed only first-order and second-order dif¬ 
ference equations. For the complete calculation 
we have encountered equations up to fourth or¬ 
der. These equations are solved with a special 
FORM program in which we substitute a suf¬ 
ficiently general combination of harmonic sums 
and then solve for the (some times thousands of) 
coefficients. The boundary values are provided 
by the Mincer program [ that was also used 
for obtaining the fixed moments. 

Problem 4 

The reduction equations have to be derived in a 
way that avoids spurious poles. Some equations 
may introduce powers of 1/e that are superflu¬ 
ous. The general rule is that one has to introduce 
at most one power of 1/e for each line that is 
eliminated. If there are more poles, the integrals 
by which they are multiplied will be needed to a 


corresponding number of extra powers of e. This 
is not always feasible. One usually avoids such 
extra poles by tedious combinations of equations. 
This is the most difficult part of the reduction 
scheme. 

Problem 5 

Next the equations have to be programmed in a 
way that produces results. Brute force applica¬ 
tion leads eventually to programs that run un¬ 
acceptably long or produce intermediate results 
that exceed the size of the available disks. The 
solution here is a very careful tabulation of inte¬ 
grals of lower complexity. For this one has to de¬ 
fine a hierarchy of complexity. Currently we have 
already more than 20000 tabulated integrals. 

Problem 6 

Because the integrals are functions of a parameter 
and involve harmonic sums of weight 6 (of which 
there are 486 different ones) each tabulated inte¬ 
gral typically takes about 20 Kbytes. This means 
that the tables are rather large, 0(500 Mbytes). 
It is not convenient to have to compile such an 
amount of code for each of the about 10000 dia¬ 
grams to be computed, even though the FORM 
compiler is very fast (about 2 Mbytes per second 
on a 1.7 GHz Pentium processor). The solution to 
this problem lies in the use of the new tablebase 
facilities of FORM [ ^ . 

Problem 7 

All topologies and subtopologies have to be pro¬ 
grammed, debugged and run. Here the old Min¬ 
cer program [ ^ for the fixed moments forms 
an indispensable tool. At any moment we can 
replace the Mellin moment A^ by a fixed inte¬ 
ger value and continue with the Mincer program. 
This way we can locate errors rather efficiently. 
Without this method the chances of obtaining the 
correct answer soon would be virtually zero. 

Problem 8 

Manage FORM [ |^. Rare bugs occur. Also 
for a problem of this complexity sometimes new 
features can bring relief. The above-mentioned 
feature of the tablebases is an example. 
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3. Status 

On the level of integrals we distinguish basic 
building blocks (BBB’s) and composite building 
blocks (CBB’s). The BBB’s are integrals in which 
the parton momentum P flows only through a 
single line in the diagram. These have been pro¬ 
grammed completely and the vast majority of 
necessary integrals have been tabulated. Occa¬ 
sionally we still have to add some integrals to the 
tables. We have reduction schemes for all CBB’s. 
They are switched on and debugged one by one 
at the moment. As there exists a rather strict 
hierarchy in the complexity of diagrams the most 
complicated will be done last, but that does not 
necessarily imply that the most complicated are 
the most difficult. The easier topologies have all 
been debugged and most of their diagrams have 
been run. At the moment the more complicated 
topologies are being treated and run. 


4. Some results 


We have completed the diagrams which con¬ 
tribute to the rij: part of the non-singlet struc¬ 
ture functions. The results for the anomalous di¬ 
mension and some parameters of the soft-gluon 
resummation are presented below. We hope to 
finish the complete non-singlet part early in the 
2003. Later that year the singlet results should 
follow. All diagrams we have run check with the 
Mincer results for several fixed moments. Our re¬ 
sults include both the anomalous dimensions and 
the coefficient functions. For the latter the reader 
is referred to ref. [ |3^ . 

The fermionic three-loop contribution to the 
even-iV non-singlet MS anomalous dimension, 
with the expansion parameter normalized as 
Q;s/( 47 r), is given by 
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Here we have suppressed the argument N of the 
harmonic sums and used the notation 

N±/(iV) = /(iV±l) 

N±i/(iV) = /(iV±*) . 

The corresponding splitting function, as usual de¬ 
fined with a relative sign, reads 
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i ln(a;)C2 - y ln(a;) ln(l+a;) + y ln^(a;) 
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where we have used 

Pqq(^) — 2(1 x) 1 X 

and all divergences for x ^ 1 are understood 
in the sense of +-distributions. The part of 
eqs. (||) and eqs. (||) has been derived before by 
Gracey [ |l| and we agree with his result. 

Our results also facilitate the determination 
of coefficients governing the soft-gluon (thresh¬ 
old) resummation [ |3^, ^ at next-to-next-to- 
leading logarithmic accuracy [ ^. 

The coefficient ^3 arising from initial-state 
collinear emissions is the coefficient of 1 /( 1 —x)+ 

/n\ 

in Pns{x). Its fermionic part reads 
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Simultaneously to our work Carola Berger [ Q 
has used a method based on eikonal expansions 
to obtain this coefficient. The calculations have 
been performed independently and put on the in¬ 
ternet before comparison. The result agrees. 

From the coefficient of ln(l—x)/(l—x)+ in the 
y part of the three-loop coefficient function we 
can furthermore determine the (complete) resum¬ 
mation coefficients B 2 and which are due 

to final-state collinear and large-angle soft gluons, 
respectively. Previously only the sum of these two 
parameters had been determined [ M from the 
two-loop coefficient function of ref. [Q. Our re¬ 
sult involves a different linear combination (with 
a prefactor f3o), and hence we can extract both 
individually, yielding 
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(5) 


The last result is intriguing and calls for further 
studies. 
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